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ABSTRACT 


This thesis addresses the critical process of assigning strike aircraft to 
targets once the targets have been identified: How do we optimally employ 
available aircraft and weapons on the current set of targets, and how can we 
modify a previously optimized assignment list to face changes in the tactical 
situation? Our contribution to the strike-planning problem includes (1) a static 
allocation model in which each aircraft makes at most one sortie during the 
planning time horizon, (2) a dynamic model in which each aircraft may make 
more than one sortie during the that horizon, and (3) extensions of these 
models with "persistence incentives,” which discourage major plan changes in 
the results when partial but important changes in the tactical situation 
necessitate reoptimization. These optimization models are mixed-integer 
programs that solve in seconds on a personal computer for realistic scenarios 
with three weapons types, 156 aircraft at seven bases, and 100 potential 
targets. In a scenario in which two new high-priority targets arise and must be 
added to an air tasking order with eight original targets, persistence incentives 
reduce the number of major plan changes from five to two. 
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EXECUTIVE SUMMARY 


The explosion of information available to decision-makers at all levels on 
the battlefield has led to tremendous strides in the ability to process, decide 
and act. At the same time, there is increasing demand for faster and better 
decisions. This thesis addresses one type of critical decision that must be made 
quickly and effectively, namely how to assign strike aircraft to targets once the 
targets have been identified. More specifically, the primary questions that this 
research attempts to answer are: 

"How do we optimally employ available aircraft and weapons on the 
current set of targets?” 

"How can we modify a previously optimized assignment list to face 
changes in the tactical situation, in a manner that balances the desire to obtain 
a new optimal solution with the desire to minimize disruption of existing plans 
based on the previous solution?” 

The general air strike planning problem is divided into five areas: target 
selection, weapon allocation, mission formation and assignment, mission 
routing and scheduling, and contingency planning. This thesis addresses 
weapon allocation, and it addresses mission allocation and assignment, for the 
aircraft of an air strike package (the "package” will typically contain other 
aircraft with non-strike roles such as suppression of enemy air defense; we do 
not directly consider these aircraft in this thesis). 

Our contributions to the strike-planning problem include (1) a static 
allocation model in which each aircraft makes at most one sortie during the 
planning time horizon, (2) a dynamic model in which each aircraft may make 
more than one sortie during that horizon, and (3) extensions of these models 
with "persistence incentives,” which discourage major changes in the results 
when partial but important changes in the tactical situation necessitate 

reoptimization. All of these models are mixed-integer programs. 

xvii 



We demonstrate these models in realistic scenarios. For instance, our 
models can allocate 156 air assets from seven bases to 100 targets in just few 
seconds on a personal computer. This particular problem has three types of 
aircraft and two possible weapons configurations for each aircraft. There are 
three types of weapons and 20 different packages can be formed from the 
aircraft available. 

The persistence paradigm plays an important role in this thesis. We show 
that a new optimal plan can differ drastically from a previous plan after only 
small changes in the tactical situation. To handle these unexpected changes 
more conveniently, the persistent model adds a term to the original objective 
function that penalizes deviations from the original plan. The new plan is 
nearly optimal in the standard sense, and unit-level planners and pilots have 
fewer plan changes to handle. 

This research is performed under the sponsorship and guidance of 
SPAWAR (Space and Naval Warfare Center). The result of this work is to be 
used inside REDS (Real-Time Execution Decision Support), a decision-support 
tool currently being developed at SPAWAR under sponsorship of the Office of 
Naval Research. The optimization modeling developed in this thesis will support 
SPAWAR's goal of having REDS provide the best solution for allocation of strike 
assets in a dynamic tactical environment. 



I. INTRODUCTION 


A. BACKGROUND AND PURPOSE 

The Space and Naval Warfare System Command (SPAWAR) in San Diego is 
developing tools for improving the planning and execution of (air) strike 
warfare. More precisely, they are working on a project called REDS (Real-Time 
Execution Decision Support), which contains two modules: the Element Level 
Planner (ELP) and the Mission Monitor (MM). ELP automates the administrative 
work involved for planning at the combat unit level and MM will monitor all 
phases of a mission, making a wealth of real-time data available for planning, 
e.g., updates of target information, aircraft load-outs, aircraft positions and 
status, etc. 

The Naval Postgraduate School’s Department of Operations Research has 
been invited to provide theoretical support in the development of the next 
stages of REDS. The effort will look for the use of the information coordinated 
by the Mission Monitor to partially automate the decision-making processes 
involved in mission planning. REDS will (i) help determine the composition of 
strike packages and the assignment of packages to targets, (ii) identify 
"efficient" (low risk, low fuel consumption, high success probability) routes for 
packages to take to targets, and (iii) provide probabilistic information on 
mobile target locations, potential actions by the adversary, etc. The mobile 
targets may be surface-to-air missiles, theater ballistic missiles, troops, etc. 

This thesis will focus on item (i), investigating different models and 
possible objective functions for optimization of strike planning. 


B. RESEARCH PROBLEM 

Since its origins in World War II, operations research methods, such as 
optimization, have dealt with the problem of resource allocation. Optimization 
has been used in many military and civilian situations to determine the best 
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assignment of tasks to agents. In these situations, optimization models ensure 
that the tasks are accomplished to the highest possible degree of performance, 
subject to constraints that must be enforced on the availability of the agents 
and the limitations of their resources. 

In our specific problem the "agents” are air strike force assets and the 
"tasks” are military targets to be neutralized. The constraints refer to 
targeting priorities, weapons effects, aircraft availability, speed, location, and 
weapon capability. The main goal is to best support the overall campaign 
objectives. Examples of this type of modeling are discussed in the literature 
review. 

We plan to go one step further than the solution of an optimization 
model for assigning agents to tasks, by developing a dynamic allocation model. 
This model allocates resources over a substantial time horizon in a changeable 
tactical environment, recommending good allocations even as real and 
contemplated modifications in the initial conditions occur. This means that the 
problems must be solved many times in rapid sequence, with the possibility of 
prior decisions constraining future ones. 

This thesis supports the strike planning effort, by partially automating 
the decision-making processes involved in mission planning. Further, it 
addresses the situation when assets are already assigned to targets and a new 
tactical scenario emerges with new targets, new priorities, aircraft 
maintenance updates, and modified weather conditions affecting weapon 
performance. 

The problem will be represented mathematically by a mixed-integer 
program (MIP). In a MIP, we have both continuous and integer variables to 
describe, quantify, and qualify the states and controls of the system. We can 
view the situation as a multiperiod problem: plans are made for multiple time 
periods in the future, and as one period elapses and better data and forecasts 
become available, the model is slid forward one period. Using ideas of 

persistence [Brown, Dell and Wood 1997], we are able to model and address 
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the very changeable tactical scenarios that arise in such real-world situations. 
And note that, even if we are planning only one period in advance, a similar 
situation can arise: A plan is developed and promulgated, but the tactical 
situation changes before that plan can be executed. Persistence is important in 
this case, too. 

C. AIR TASKING ORDERS AND AIR STRIKE PACKAGES 

An Air Tasking Order (ATO) is the administrative vehicle used to 
disseminate daily plans to units and to command and control agencies. The ATO 
normally provides specific instructions to include radio call signs, times on 
targets, and other detailed information required for the execution of a plan. An 
ATO is the result of a complex process of target selection and allocation of 
assets covering a myriad of missions. It usually takes two days of planning and 
strict control of execution to develop an ATO. ATOs require careful 
coordination of tasks and consideration of weaponeering data, force structure, 
sortie and weapons availability, intelligence aspects of the enemy, weather, 
and numerous other pieces of information. 

There are usually three ATOs in existence at any one time: (1) the ATO 
in execution (today’s plan), (2) the ATO in production (tomorrow’s plan), and 
(3) the ATO in planning (following day’s plan). The idea of reducing that cycle 
is represented in Figure 1 according to the concept of a real-time C 2 
infrastructure. 

Griggs [1994] precisely defines an air strike package: "An strike package 
is a group of fighter and bomber aircraft that have combined to provide mutual 
support against enemy threats while they achieve a common goal of destroying 
a set of targets. Strike packages are normally constructed in several steps. 
First, the mission planner must select the right type and number of aircraft and 
munitions to efficiently destroy each target. Next, all flights attacking targets 
in the same vicinity are grouped into packages if aircraft speed restrictions and 
tactics are compatible. Last, the mission planner must add suppression of 
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enemy air defense (SEAD) aircraft and air-to-air fighter escort, or sweep, 
aircraft to protect the groups of attackers.” 

The addition of SEAD and escort aircraft depends on their availability, 
the enroute threats such as surface-to-air missiles, the mission, and the type of 
aircraft in the package. A group of flights attacking targets in the same vicinity 
together with SEAD and escort aircraft comprise a typical air strike package. 
Figure 2 shows an example of a strike package used during the Gulf War. 


Detect 


Decide 


Deliver 


Kill 


BDA 


Day 1 

Day 2 

Day 3 

Day 4 

Day 5 

Day 6 



Figure 1. ATO timelinei 


The actual production cycle of an ATO (Air Tasking Order) 
comprises three days of preparation and execution, but the 
increasing flow of information from the battlefield is driving the 
need for a more compressed timeline. 


1 The source of this figure is the Naval Aviation Interoperability Assurance Office, 
(http://www.dtic.mil/ndia/systems/Clark.pdf) "Streamlining Acquisition Through Collaborative 
Engineering,” 10/06/2000. 
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Strike force planning can quickly become very difficult to model and 
solve if we consider the huge number of possible combinations of aircraft, 
tactics, and weapons against each type of target, threat and environmental 
condition. In order to derive the best methodology for our specific objectives, 
we first review the current status of research on strike force allocation. 



Figure 2. Gorilla Package2 

A "gorilla package” was a type of strike package that 
placed a large number of aircraft over a target in a short period 
of time during the Gulf War. Planning time and effort increased 
dramatically with the size of the package. 


2 Cohen 1993, Vol IV, pp. 164-165. 
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There are different problem characteristics that may be modeled 
depending on the model’s purpose. The allocation of strike forces can be either 
part of a real-time planning tool or a long-term budget planning system. Each 
problem focuses on different aspects of mission formation and assignment. 
Acceptable solution time depends on the model’s purpose. 


D. LITERATURE REVIEW 

Our research sponsor at SPAWAR, John McDonnell, introduced us to a 
model for composing strike packages from available assets and allocating those 
packages to targets. He and colleagues Gizzi and Louis [McDonnell et al. 2001] 
develop a nonlinear optimization formulation of the allocation problem, and 
solve it approximately using a genetic algorithm (heuristic). Their approach 
encompasses both strike and suppression responsibilities as well as multi-target 
and multi-threat allocations. The model assumes that weaponry load-outs have 
been predetermined and the strike force can be reassigned to new targeting 
objectives. 

Abrahams [1998] and Balart [1996] are the most important reference for 
the work of McDonnell et al. Abrahams and Balart’s model is based on a non¬ 
linear objective function for a static allocation, where crucial parameters such 
as target value and weapon effectiveness are independent of time. A genetic 
algorithm provides approximate solutions. 

Li, Curry and Boyd [2002] are currently working on an integer 
programming model for the strike force asset allocation problem, sponsored by 
the Office of Naval Research. Although this report focuses on determining if a 
solver can handle a large integer model in real time, it contains interesting 
insights on how to build suppression packages to act against threats. 

Criteria for building strike packages are found in the USAF models 
described by Griggs [1994], Jackson [1989], and Yost [1995]; although each 
author uses a different objective and different interpretations for "sortie,” 
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"target,” "weapon,” etc. Our research incorporates some of the main ideas 
from these Air Force models. 

According to Yost [1996], the Time Strike munitions optimization model 
was introduced in 1995 for use by various USAF agencies to develop 
requirements for conventional munitions, to refine operational plans based on 
the availability of different mixes of munitions, and to assess the effects of 
procuring different types and quantities of munitions. Time Strike’s objective 
pursues the phases (sets of goals for each target class) in a hierarchical order 
defined by the user. The model’s notion of target classes is a major difference 
from previous models, and supports the fact that campaign objectives involve 
killing collections of related targets rather than individual targets. It is a multi¬ 
period linear programming model. 

Koewler [1999] creates a prototype scheduling and allocation tool that 
strikes a balance between ease of use, accurately defining and solving the 
allocation problem, and generating solutions in an operationally acceptable 
amount of time. A combination of concepts from project scheduling, object- 
oriented programming and heuristics form the basis for the methodology 
developed. 

Saling [1999] examines the Joint Air Operations Center structure, the 
relationship of information to the Master Air Attack Plan and methods of 
distributing that information to the warfighter through the Air Tasking Order 
and alternately through dynamic re-tasking. 

Dolan [1993] offers a solution to the problem of producing a timely and 
flyable ATO (Air Tasking Order) that effectively uses assigned aircraft. His 
model decides which strike package should be assigned against each target and 
which available launch sites should provide the assets required in the selected 
strike packages. The output is an ATO used in theater level war games 
conducted at the Naval War College. 
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Crawford [1994] enhances a version of Dolan’s optimization model to 
explicitly incorporate the dimension of time, thereby allowing multiple sorties 
per aircraft per day, something not allowed in Dolan's model. Crawford’s model 
takes packages as inputs, which are provided by a team in a war game 
exercise. 

The doctrine currently used by American military forces completely 
guides the process of mission planning. Joint Publication 3-01.2 [1986] provides 
us some insights into planning operations, discussing the Intelligence activity of 
targeting, as well procedures and techniques of Command and Control. This 
document describes some fundamental principles from which we define our 
objectives. Two of them are especially important when conducting counterair 
operations: the concentration of force and the economy of effort. 

Concentration of force is the effective application of combat power, 
which requires that sufficient force be concentrated at the appropriate time 
and place to ensure achievement of the objective. 

Economy of effort is the correct selection and use of weapon systems, 
sound distribution of forces and careful balance in the allocation of tasks. 
When applying this principle, the commander intends to achieve effective 
concentration of power at the decisive time and place while conserving 
weapons for countering enemy reattacks. 

The same document also instructs on how to establish target priorities. 
Five criteria are of vital importance: threat, feasible effect (degree of positive 
effect, in terms of degrading every capability or enhancing friendly 
operations), delay in effect (time between the initial engagement and the 
desired effect; concentration of effort may compress that time), risk 
calculation, and forces available. 

Part of our research is directed to answer the question of how to modify 
a previously optimized plan. We are particularly interested in solving a 
multiperiod model in a rolling-horizon format, where the model is solved at the 
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beginning of some period or by suggestion of the decision maker. The 
recommendations of the solution for the previous period are remembered and 
deviations from them are charged a penalty in the objective function. This is 
called a "persistence incentive.” 

The method we use for handling persistence derives from Brown, Dell 
and Wood [1997]. These authors state that the lack of persistence is one of the 
most common sources of complaints about optimization. They address the issue 
using a series of case studies that demonstrates how persistence can mediate 
the differences in focus between managers and modelers, and show how to 
develop models from the start with persistence in mind. 


E. PHASES OF AIR STRIKE PLANNING 

The general air strike planning problem can be structured to identify and 
classify objectives, constraints, decisions, and influencing factors. Using this 
structure, decisions are divided into five problem areas: (i) target selection, 
(ii) weapon allocation, (iii) mission formation and assignment, (iv) mission 
routing and scheduling, and (v) contingency planning [Glenn 1980]. 

The target selection activity examines potential targets to determine 
military importance, priority of attack, and weapon feasibility to obtain a 
desired effect [USAF 1998]. The selected target systems are then further 
analyzed to determine their components and critical elements. This phase 
distills the commander’s objectives into a list of targets. The product of this 
phase is a suggested target list with recommended priorities assigned and the 
extent of desired damage specified. 

Weapon allocation (also called "weaponeering”) estimates the quantity 
of a specific weapon type required to achieve the desired level of damage to a 
given target, considering target vulnerability, weapon effects, munitions 
delivery errors, damage criteria, weapon reliability, etc. [USAF 1998]. Weapon 
effectiveness varies according to the weapon, target, damage criteria, delivery 


9 



conditions, and target environment. There are different ways of stating 
weapon effectiveness according to the target/weapon combination. We will 
derive and use only the probability of kill, based upon information about the 
target, the weapon and the aircraft available. 

Mission formation and assignment designs the actual strike package. 
Usually, planners start with the allocation of strike assets and then assign the 
non-attack mission platforms (i.e., aircraft to escort, suppression of enemy air 
defenses, jamming, airborne control, tactical reconnaissance, air refueling, 
and search and rescue), which support the package’s ingress and egress. 

During the mission routing and scheduling process, planners ensure that 
the mission packages are "deconflicted” with other mission packages, tankers 
are available at the refueling points, times of launch and landing are 
synchronized, etc. The deconfliction process prevents conflicts such as aircraft 
on different missions using the same altitude or too many aircraft requiring 
refueling or landing at the same time. This coordination must still conform to 
precise arrivals at the required time on target. The decision-maker also 
considers the enemy’s air defense and weather conditions enroute in this 
phase. 

Finally, contingency plans specify secondary targets and conditions for 
diverting the strike to them, and variations in strike tactics if the weather 
conditions at the target differ from those anticipated. 

This thesis addresses two areas of strike mission planning: (ii) weapon 
allocation and (iii) mission allocation and assignment. Only the strike 
components of the packages are considered. Additional work is required to 
include specialized escort and SEAD aircraft to increase the package’s 
probability of success against the enemy’s defenses. 
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1. Weapon Allocation 

In the weapon allocation problem, we must consider the effectiveness of 
a specified number of weapons delivered by a given number of aircraft against 
a target that has an associated minimum damage criterion or minimum 
acceptable probability of kill. To accurately model this situation requires 
nonlinear functions for evaluation of the probability of killing the target. We 
would prefer, however, a linear optimization model for assigning strike 
packages to targets (linear in the optimization model’s decision variables and 
constraints). So, how do we incorporate a nonlinear measure of effectiveness in 
a linear optimization model? 

The solution is to enumerate the "strike options,” i.e., possible strike 
packages, and calculate their effectiveness against each target as part of the 
mixed-integer programming (MIP) model’s input, rather than as part of the 
model’s decision process. Then the MIP can be defined with binary variables 
representing the possible pairings of targets to strike packages. Dolan’s [1993] 
and Crawford’s [1994] models also incorporate pre-enumerated strike 
packages, with a subjective preference rating for each given target type. 
Ratings from one to five are specified, with one being the most preferred and 
five the least. 

Instead of a preference rating, we calculate the probability of kill for 
each package against each target type. This calculation is simplified and only 
depends on the type and quantity of weapons the aircraft in the package can 
launch. We also consider that some weapons don’t work against certain types 
of targets or under certain weather conditions. 

Each possible combination of weapons load for a given aircraft is called 
a configuration. An example of a strike package could be the assignment of 
two aircraft of type a in configuration c, where configuration c contains four 
weapons of type w 1 and two weapons of type w 2 . 
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2. Mission Allocation and Assignment 

The phase of strike planning that receives the most emphasis in this 
thesis is mission allocation and assignment. It is extremely difficult to address 
the problem with full generality. Therefore, choices must be made concerning 
alternate ways of handling various aspects of the problem. As a result, there 
could be many reasonable modeling approaches. 

The first choice to be considered is whether the model is static or 
dynamic. A static model treats all strikes as taking place during a single time 
period. This does not mean the strikes are literally simultaneous, but the 
resolution of modeling ignores time effects. Thus, the static model does not 
explicitly consider the possibility of aircraft completing a strike, returning to 
its launch site, refueling, and launching for a second strike. The second launch 
would have to be considered in a second run of the static model. 

A dynamic model, on the other hand, explicitly treats the passage of 
time, allowing aircraft to perform multiple strikes. The dynamic model can also 
include the feature of having the model decide when targets should be struck 
within a given time window. This thesis develops both static and dynamic 
models. 

A second modeling choice is whether sorties must strike only one target 
or are allowed multiple targets. The single-target restriction may not be 
entirely realistic, but is much easier to model than the multiple target case. 
Throughout the thesis, the single-target restriction is enforced. 


Table 1. Summary of modeling options 


Time 

Static vs. Dynamic 

Targets per sortie 

Single target vs. Multiple target 

Strike package composition 

Homogeneous vs. Heterogeneous 
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Finally, strike packages can be homogeneous or heterogeneous with 
respect to the types of strike aircraft in the package and with respect to the 
types of weapons carried in each aircraft. The static model developed in this 
thesis allows for heterogeneous strike packages. For computational reasons, 
the dynamic model is restricted to strike packages consisting of only one type 
of strike aircraft, each with the same weapons loadout. Table 1 presents a 
summary of modeling options for the air strike asset allocation problem. 

3. Contingency 

Although the principal emphasis of a strike plan must be on events that 
are judged most likely to occur, it is also important for the plan to offer 
responses to the less likely possibilities [Glenn 1980]. A strike plan should 
include contingency plans that provide secondary targets, backup assignments 
and variations in strike tactics if weather conditions, aircraft availability or 
target characteristics change. 

Our models do not directly provide contingency plans, in contrast to 
Kuykendall’s [1998] Tomahawk missile assignment problem, which assigns 
backups for every mission. However, our models can be helpful for dealing with 
contingencies indirectly. As conditions change, the models can be re-run. The 
persistence feature will encourage the new solution to be similar to the 
previous one, making it easier to adapt to the necessary changes. 

The models described in this thesis incorporate important aspects of the 
strike planning process to produce air tasking orders. Chapter II describes a 
static model. Chapter III treats a dynamic case. Chapter IV develops extensions 
of the static and dynamic models with persistence incentives. Finally, 
Chapter V summarizes our results and suggests directions for future research. 
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II. STATIC MODEL WITH HETEROGENEOUS PACKAGES 


This chapter develops a static optimization model for creating air 
tasking orders. The key decisions are: which combination of aircraft, with 
which weapon loadouts, coming from which launch site, should strike which 
target? 

The mixed-integer program makes the best overall assignment of strike 
packages, taking into account the most important characteristics of the assets 
and targets. Key features of this model are given in Table 2. 

The model in the next chapter is dynamic. Although some parameters 
and variables are similar for both models, the models are defined completely 
and independently in separate chapters. 


Table 2. Static model summary 


General 

A MIP (mixed-integer program) with a static time 

approach allocates heterogeneous strike packages, 

allowing only one target struck per sortie 

Objective Function 

A linear combination of three objectives: minimize 

value of targets not assigned, minimize the effects of 

imperfect matching of targets to packages 

(incomplete damage, long-distance flight, etc.), and 

maximize value of unused aircraft 

Campaign Objectives 

Input as a required time on target (TOT) 

Target Prioritization 

Based on target values 

Aircraft Packaging 

Different types of aircraft carrying different mixes of 

weapons are allowed (not restricted to homogeneous 

packages) 
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A. INDICES 

The indices used to define this model are: 


a 

aircraft type 

{’FA18’, ’A-6’, ’A10’,...} 

c 

configuration 

{’C 1’, ’C2’, ’C3’,...} 

i 

launch site 

{’Airbase-01’, ’CV71’,...} 

j 

target 

{’Safwan’, ’SCUD’,...} 

k 

target type 

{’Airfield’, ’Bridge’,...} 

n 

strike package 

{’N1 ’, ’N2’,...} 

r 

threat type 

{’AAA’, ’SAM’, ’Air-to-Air’,...} 

w 

weapon 

{’MK84’, ’GBU31 ’, ’AGM65’,...} 


Launch sites can be either airbases or aircraft carriers. With this 
understanding, the terms "site” and "base” are used interchangeably. An 
aircraft configuration is an instance of a possible loading of guns, ammunition, 
missiles, bombs, and external fuel tanks on a single aircraft, although this 
thesis ignores fuel tanks for simplicity. Potential configurations are specified in 
a descriptive document called standard conventional loads (SCLs). An example 
for the F/A-18 aircraft is given in Table 3: 


Table 3. Example of configuration 


Config ID 

MG25 

Type A/C 

F/A-18 

Gun 

20mm 

500 rounds 

Missile 

AIM7 

2 

Missile 

AIM9 

2 

Bomb 

MK82 

5 

External fuel 

TANK 

2 
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The strike package is a mixture of different aircraft and configurations. 
Table 4 gives an example: 

Table 4. Examples of heterogeneous strike packages for the static model 



The strike packages are assumed to be given information for the MIP. In 
practice, they are obtained from past experience and/or software that 
enumerates feasible combinations and evaluates them with weapons 
effectiveness models such as in appendix A. 

In attacking a target, strike aircraft will typically be subject to threats 
from surface-to-air missiles (SAM), anti-aircraft artillery (AAA), etc. Mission 
success will depend on surviving these threats so they are a part of our models. 
However, we do not explicitly model the auxiliary aircraft in a strike package 
that might be responsible for dealing with these threats. 

B. PARAMETERS 

The data parameters used to define this model are: 

1. Asset Data 

acval(a) value of aircraft a if preserved for later use 

acavail(/,a) number of aircraft a available at site ; 

range(a,c) range in nautical miles of aircraft a in configuration c 

psv(a,c,r) probability of survival against threat type r for 

aircraft a in configuration c. This will be increased 
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when considering the use of support missions (escort 
and SEAD) in the same package 


2. Target Data 

tthreat (J,r) 
mindamage(;) 

tgtvalO') 


equals one if threat type r is present at target j 
minimum required damage level on target j 
target value 


The minimum required damage on a target is sometimes specified with 
subjective terms like "light”, "moderate” and heavy.” Throughout the thesis, 
the minimum required damage on a target is interpreted numerically as a 
specified probability of kill (pk). The mapping between the subjective 
terminology and probabilities is given in Table 5. 


Table 5. Levels of damage and probability of kill 3 


Damage Level 

Description 

pk 

Light 

Minor damage, some functions lost, but 
still capable of operation 

0.3 

Moderate 

Extensive damage, many functions lost. 
Operation still possible but at reduced 
effectiveness 

0.7 

Heavy 

Unable to operate 

0.9 


3. Strike Package Data 

numweap(a, c,w) number of weapon w carried by aircraft a in 
configuration c 

numac(n,a,c) number of aircraft a in configuration c flown in 

3 Class notes from Introduction to Naval Weapons Engineering (ES 310), Damage Prediction 
(http://www.fas.org/man/dod-101 /navy/docs/es310/dam_crit/dam_crit.htm) 
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pk (n,k) 


strike package n 

probability of kill for package n against target type 
k. The method for computing pk prior to the 
optimization is given in Appendix A 

refuelmax maximum number of air refuelings per strike 

minpk minimum probability of kill. Used to ensure that a 

package without sufficient damage capability will 
not be assigned to any target 

4. Geographic Data 

dist(ij) distance in nautical miles between site ; and target j 


5. Derived Data 

nac(n,a) number of aircraft a in package n 

nw(n,w) number of weapon w in package n 

refuels(n,/J) number of refuelings needed by package n from 
site ; when attacking target j. The aircraft in the 
package with the lowest range will dictate the 
value of this parameter 


nac (n,a) = ^]numac(a,c) 

c 

nw(n, w) = ^numac(n,a,c) ■ numweap(a, c,w) 

a,c 


refuels(n,/,;') = max 

(a,c) s.t. numac(n,a,c)>0 


| dist(?,j) j 
[|_range(a,c)Jj 


dpen (n,i,j) distance penalty when package n from site ; strikes 

target j. Following Crawford [1994], the distance 
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penalty increases as the length of the mission 
approaches the range of the aircraft. When a combat 
radius is exceeded, the aircraft must refuel which 
causes a jump in the distance penalty 

ml proportionality constant of dpen for distance 

m2 proportionality constant of dpen for refueling 

dpen(n,/,;) = ml ■ max 

(a,c) s.t. numac(n,a,c)>0 

+ m2 ■ refuels(n,/,;') 

pkpen(nj) penalty imposed if package n’s probability of kill is 
different from the required damage for target j 

m3 proportionality constant of pkpen when the probability 

of kill is not enough 

m4 proportionality constant of pkpen when the probability 

of kill is greater than that required 


dist(/,;) 

range(a,c) y 


pkpen(n, j) = m3 ■ 


mindamage(;) - pk(n, j) 
mindamage(;) 


pkpen(n, j) = m4 ■ 


pk(n, j) - mindamage(;) 
mindamage(;) 


if mindamage(;') > pk(n,;) 
if pk(n,;) > mindamage(;) 


attrition(nj) expected attrition if package n strikes target j 
m5 proportionality constant of attrition 

attrition^,;) = m5 ■ ^ numac(n,a,c) ■ 1- ]^[ psv(a,c,r) 

(o,c) s.t. numac(n,o,c)>0 _ retthreat( 7 ,r) 
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6. Objective Function Data 

stress (n,i,j) penalty for imperfect damage, long distance 
(refueling) and low probability of survival against 
enemy defenses if package n is flown from site ; 
against target j 

stress (n, i, j) = dpen(n, i, j) + pkpen(n, j) + attrition(n, j) 


C. DECISION VARIABLES 

The primary decision variables are binary. They allow for selection of 
which package is assigned to each target and which site provides the assets 
required. 


fl 


STRIKE (n,/,;) = 


10 


if strike package n is assigned from 

site / to target j 

otherwise 


This variable has its domain restricted to implicitly enforce constraints 
and to make solution easier. In particular, STRIKE(n,/J) exists only if all the 
following conditions hold for (n,i,j): 

refuels(n,/,;) < refuelmax 

nac(n,a) < acavail(/,a) Va s.t. nac(n,a)>0 

pk(n,/,;') > minpk 


There are two additional sets of variables: 

NOGO(;) = 1 if target j not attacked, zero otherwise 

PRESERVER, a) number of aircraft from site ; not assigned to any target 

These variables have integer interpretations but can be treated as 
continuous variables: by virtue of the constraints’ structure, they cannot 
fractionate in any feasible solution. 
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D. OBJECTIVE FUNCTION 

The objective function has three components which may be weighted or 
solved for lexicographically: 

Minimize ^tgtval(y) • NOGO(;) 

] 

Minimize ^stress (n,i,j) ■ STRIKE(n,/,;) 

n,i,j 

Maximize ^acval(a) ■ PRESERVER, a) 

/, a 

The first component term charges a penalty when a target is not struck. 
The second component accounts for imperfect damage, long distance 
(refueling), and lower probability of survival against enemy defenses. The last 
one rewards savings of aircraft for use in future missions and unforeseen 
contingencies. 

Throughout this thesis, multiple objective-function components are 
combined into a weighted sum, with positive weights on the terms to be 
minimized and negative weight on the terms to be maximized. It is also 
possible to do a lexicographic optimization in which the highest priority 
objective is optimized first; then it is constrained to its optimal value while the 
second priority objective is optimized; and the process continues as long as 
there are alternative optima at each stage [Rosenthal 1985]. 


E. CONSTRAINTS 

The static model needs only two sets of constraints. The first set ensures 
that the targets are struck or penalties are incurred through the variables 
NOGOtf): 

£STRIKE(n,/,;) + NOGO(;) = 1 V; 

n,i 

The second set of constraints ensures that we only use aircraft that are 
available: 
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^nac (n,a) ■ STRIKE(n,/,;) + PRESERVER, a) = acavail(/,a) Vi,a 

n,j 


F. RESULTS 

We have implemented the static model using GAMS with CPLEX and XA as 
solvers. The results for a 100-target problem are given in Table 6. 


Table 6. Static model results for 100 targets 


Running on a Dell Workstation Precision 340 (Pentium IV 2 GHZ 1 GB RAM). 
GAMS IDE environment version 2.08.3, revision module 117, October 16, 2000. 



CPLEX 

(version 6.6.1) 

XA 

presolver 

1 row and 9 columns 
eliminated 

- 

problem size 

108 rows, 5900 columns, 
12300 nonzeros 

109 rows, 5909 columns, 

18209 nonzeros 

optcr = 0.05 

< 1 sec 

< 1 sec 

optcr = 0.0 

< 2 sec 

< 3 sec 


This problem also has seven bases with three types of aircraft and two 
possible configurations for each aircraft. There are three types of weapons and 
20 different packages are allowed to be formed from the total of 156 aircraft 
available. 

Observe that CPLEX applies a "presolver” phase, which reduces the size 
of the MIP. The parameter "optcr” is a relative measure of optimality, a bound 
on how far from the best possible answer we are; "optcr = 0.05” means we are 
no more than 5% off. The smaller the optcr, the more time is needed for the 
solver to find a solution. 
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III. DYNAMIC MODEL WITH HOMOGENEOUS PACKAGES 


This chapter presents the dynamic extension of the static optimization 
model in the previous chapter. Adding the time dimension allows explicit 
consideration of assigning aircraft to multiple sorties. In exchange for this 
added realism, the dynamic model introduces the restriction of homogeneous 
strike packages. Each package contains only one type of strike aircraft, and 
each aircraft in a package carries the same type of weapon. The model is 
summarized as follows (Table 7): 


Table 7. Dynamic model summary 


General 

A MIP allocates homogeneous strike packages to 

targets over a multi-period time horizon 

Objective Function 

Linear combination or lexicographic optimization of 

four objectives: minimize target value of targets not 

assigned, minimize attrition, minimize distance 

penalty, and maximize value of unused aircraft 

Campaign 

Objectives 

Input as a required time on target and minimum 

damage for each target. Target prioritization based on 

target values 

Aircraft Packaging 

Strike packages are created for valid combinations of 

aircraft and configurations, but are restricted to be 

homogeneous. Packages are formed by aircraft of the 

same type, from the same site, with the same 

configuration. We also require that each aircraft 

deliver the same amount of only one type of weapon 

Planning Horizon 

The user specifies the present time, time of the first 

period, time between periods, and total number of 

periods 
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A. INDICES 


a 

aircraft 

{'FA18’, f A-6’, 'A10’,...} 

c 

configuration 

{'Cl’, f C2’, f C3’,...} 

i 

site 

{ f Airbase-01 ’,...} 

j 

target 

{ f Safwan’, f SCUD’,...} 

k 

target type 

{'Airfield’, 'Bridge’,...} 

n 

strike package 

{'N1 ’, 'N2’,...} 

nac 

number of aircraft 

{'NAC1 ’, 'NAC2’,...} 

nw 

number of weapons 

{'NW1’, 'NW2’,...} 

r 

threat type 

{'AAA’, 'SAM’, 'Air-to-Air’,...} 

t 

time period 

{'T1 ’, 'T2’, 'T3’,...} 

w 

weapon 

{'MK84’, 'GBU31 ’, 'AGM65’,...} 


The index nac is a device for representing general integer variables as 
binary variables. For example, if a binary variable for a package with nac = 
"nac4” is positive, then there are four aircraft in the selected package. 
Although contemporary integer programming solvers handle general integer 
variables, we use this device to allow for the calculation of the nonlinear 
probability of kill prior to optimization. The index nw is similar to nac but 
counts weapons per aircraft, rather than aircraft per package. 


B. PARAMETERS 

1. Asset Data 

numweap(a,c,w) number of weapons of type w carried by aircraft a 
in its configuration c 

cep(a,w) circular error probable (in feet) of weapon w 
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lethal(w,k) 


wstate(w) 

psv(a,c,r) 

speed(a,c) 

acval(a) 

recover(a,c) 

range(a,c) 

acavail(/,a) 

dueback(/,a,t) 


delivered by aircraft a. This is interpreted as the 
radius of a circle within which half of the delivered 
weapons are expected to fall 

lethal radius (in feet) of weapon w against target 
type k. Beyond the lethal radius, we consider the 
weapon to have no effect, whereas within the 
lethal radius it has constant effect. This "cookie- 
cutter” approach is the conceptually simplest kind 
of weapon effect model 

minimum weather state in which weapon type w 
can operate. Equals one for "poor,” two for 
"marginal” and three for "good” 

probability of survival against threat type r for 
aircraft a in configuration c. This will be increased 
when considering the use of support missions 
(escort and SEAD) in the same package 

speed (in knots) of aircraft a in configuration c 

aircraft value 

recovery time (in minutes) of aircraft a to 
configuration c, includes all services (reloading, 
refueling, etc.) needed to prepare the aircraft for 
the next sortie upon return from another 

range (in nautical miles) of aircraft a in 
configuration c 

number of aircraft of type a available at site ; 

number of aircraft of type a from site ; to be 
available in time period t. For the first time period, 
this is the number of planes initially available. For 
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dist(/J) 

pack (n,a,c,nac) 

refuelmax 

minpk 

2. Target Data 

tthreat (J,r) 
mindamageO') 

tgtvalO') 

reqtotO') 
forecast (t) 


later time periods, it refers to aircraft from a 
previous plan that have not yet returned from 
earlier sorties 

distance (in nautical miles) from site ; to target j 
(default is "great circle” if other data are not 
available) 

number of aircraft a with configuration c in 
package n. The packages are generated considering 
the realistic combinations of a, c and number of 
aircraft 

maximum number of air refuelings for all packages 

minimum probability of kill. Used to ensure that a 
package without sufficient damage capability will 
not be assigned 


equals one if threat type r is present at target j 

desired damage level on target j, i.e., the minimum 
acceptable probability of kill for the strike package 
to be assigned to target j 

target value: this information will be important to 
derive objective function coefficients 

required time on target (TOT) 

weather state forecast for period t. Equals one for 
"poor,” two for "marginal” and three for "good” 
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Figure 3. Time related parameters 


Time in the dynamic model is discrete, which means we 
divide the planning cycle into a finite number of "periods.” The 
planning of an ATO occurs at present time "to.” Considering that 
some time is needed to process and disseminate the orders, we 
define a "gap” before the first time "tj” when a launch can 
occur. We also define "resolution” as the interval between 
successive periods and "horizon” as the total cycle length, from t ? 
to the end of the last period. 


3. Time Control Parameters 

The time related parameters described in Figure 3 are the following: 
pt present time (t 0 ) 

gap gap interval (t ? - t 0 ), the minimum time required for 

an order to be completely processed 

hz horizon (t/ as t), the number of time periods 

rs resolution (t r - t 2 ), length of each time period 


4. Derived Data 

pk (n,k,w,nw) probability of kill for package n with nw weapons w 
from each aircraft against target type k. Refer to 
appendix A for the calculations of this parameter 

trvtime(/, a,c,j) number of periods that aircraft a in configuration c 
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will take to go from site / to target j 

fit (i,a,c,j) time (in minutes) required for a round-trip mission 

including recovery of aircraft a in configuration c 
from site ; when attacking target j 

flp(;,a,cj) number of periods required for a round-trip mission 
including recovery of aircraft a in configuration c 
from site ; when attacking target j 

refuels(/,a,cj) number of refuelings needed by aircraft type a with 
configuration c from site ; when attacking target j 


trvtime(/,a,c,;') = 


60 ■ dist(/, j) 
rs ■ speed(a,c) 


flt(;, a, c, j) = 2 ■ 60 dist(/,y) + recover ( Q) 
speed(a,c) 


f[p(i,a,c,j) = 


fit (i,a,c,j) 


rs 


refuels(/,a,c,;') 


dist (i,j) 
range(a,c) 


dpen(/,a,cj) distance penalty when aircraft a in configuration c 
from site i strikes target j 

ml proportionality constant of dpen for distance 

m2 proportionality constant of dpen for refueling 

dpen(/, a, c, j) = ml ■ — + m2 ■ refuels(/, a, c, j) 

rang e(a,c) 

pkpen {n,j,w,nw) penalty imposed if package n’s probability of kill is 
different from the required pk for target j 

m3 proportionality constant for pkpen when the 
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m4 


probability of kill is less than required 

proportionality constant for pkpen when the 
probability of kill is greater than that required 


if mindamage(7) > pk (n,j,w,nw ): 

, , . . mindamagei/') -pk(n, j, w.nw) 

pkpen(n, j, w, nw) = m3-——— K v — — -- 

mindamage(;) 

if pk(n,;, w,nw) > mindamage(7) 

pk(n,7,w,nw) - mindamage(7) 


pkpen(n, j, w, nw) = m4 


mindamage(;) 


dnacc(n,a,c) 

dnac(n,a) 


number of aircraft a in configuration c flying 
with strike package n 

number of aircraft a in package n 
dnacc(n,a,c) = ^nac ■ pack (n,a,c,nac) 

nac 

dna c(n,a) = ^dnacc(n,a,c) 


attrition (n, j) 


expected attrition if package n strikes target j 


attrition(n, j) = nac 


1- Y[ psv(a,c,r) 

re tthreat (j,r) 


\/a,c,j,nac 


pot(;,t) equals one if required time on target for target j 

is during period t, zero otherwise. Derive pot 
according to: 

pot(;, t) = 1 if rs • (t -1) < reqtot(;) - pt - gap < rs ■ t 
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c. 


DECISION VARIABLES 


M 


STRI KE(n, i, j, w, nw, t) = 


[0 


if package n coming from site i is assigned 
to strike target j with nw weapons of type w 
from each airplane with take-off during period t 
otherwise 


The binary variable STRIKE (n,i,j,w,nw,t) exists only if the following 
conditions hold: 

pot(;,£ + trvtime(/,a,c,;')) = 1 Va,c s.t. dnacc(n,a,c) > 0 
refuels(/,a,c,;) < refuelmax Va,c s.t. dnacc(n,a,c) > 0 
assign(/,a) > 0 Va s.t. dnac(n,a) > 0 

wstate(w) < forecast(f) 
pk(n,;',w,nw) > minpk 


NOGO(;) equals one if target j is not attacked, zero 

otherwise (continuous variable) 

PRESERVE(/,a,£) number of aircraft a from site ; not assigned to 

any target in period t 

The last two variables have integer interpretations but can be treated as 
continuous. By virtue of the constraints’ structure, the variables cannot be 
fractional in any feasible solution. 


D. OBJECTIVE FUNCTIONS 

min ^tval(;') • NOGO(;') 

j 

max ]T acval(a) ■ PRESERVER, a, t) 

i,a,t 

Z [dpen(i,a, /')+pkpen(n, w,nw)+attrition(n, /')1 ■ 

n,i,a,j,w,nw,te dnac(n,a) ' STRIKER, 7 , j, W, HW, t) 
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E. CONSTRAINTS 

Strike all targets or don’t 

X STRI KE(n, i,j, w, nw, t) + NOGO(;) = 1 V; 

n,i,w,nw,t 


Balance of aircraft 

^ dnac(n,a) ■ STRIKE(n,/,;,w,nw,t) + PRESERVE (i,a,t) = 

n,j,w,nw 

dueback(/, a, t) + PRESERVER, a, t -1) + 

^ dnacc(n,a,c) ■ STRIKE(n,/,;', w,nw,t - flp(z, a,c,j)) V i,a,t 

n,c,j,w,nw 

The left-hand side of a balance constraint is the number of aircraft of 
type o that launch from site / during period t or remain on the ground for 
future use. The right-hand side of the balance constraint is the sum of 
previously preserved aircraft and those returning from earlier mission during 
period t. Earlier missions represented by the dueback parameter were 
launched prior to t h whereas the longer term on the right-hand side represents 
returning missions launched after t ? . 


F. RESULTS 

We have implemented the dynamic model using GAMS with CPLEX and XA 
as solvers. The results for 100 targets and 16 periods of 15 minutes each are 
given in Table 8. 

This problem has seven bases with three types of aircraft and two 
possible configurations for each aircraft. There are three types of weapons and 
24 different packages are allowed to be formed from a total of 36 available 
aircraft. 


33 



Table 8. Dynamic model results for 100 targets 


Running on a Dell Workstation Precision 340 (Pentium IV 2 GHZ 1 GB RAM). 
GAMS IDE environment version 2.08.3, revision module 117, October 16, 2000. 



CPLEX 

(version 6.6.1) 

XA 

presolver 

1213 rows and 1219 columns 
eliminated 

- 

problem size 

98 rows, 10093 columns, 
20094 nonzeros 

1361 rows, 11362 columns, 
32661 non-zeros 

optcr = 0.05 

< 1 sec 

< 1 sec 

optcr = 0.0 

201 sec 

> 1200 sec 
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IV. ALLOCATION WITH PERSISTENCE 


Mathematical programming models can be stated and solved so that they 
exhibit varying degrees of persistence with respect to previous values of 
variables, constraints, or even exogenous considerations [Brown, Dell and Wood 
1997]. The importance of the persistence emerges when small changes to 
input data lead to drastically different solutions. 


A. PERSISTENCE FOR THE STATIC MODEL 

The following scenario provides an example of how to achieve better 
results by considering persistence. Initially, we have a list of ten targets to be 
struck and a limited number of aircraft available. Table 9 shows the possible 
combinations of packages: 


Table 9. Packages for the static model’s scenario 


package 

aircraft 

configuration 

number of a/c 

nOI 

A-10 

cl 

2 

n02 

A-10 

cl 

4 

n03 

A-10 

cl 

2 

A-10 

c2 

2 

n04 

F-16 

cl 

1 

F/A-18 

cl 

1 

n05 

F16 

c2 

1 

F/A-18 

c2 

1 

n06 

F16 

cl 

2 

F/A-18 

cl 

2 


The configurations differ with respect to weapon loading, as described in 
Chapter II. These aircraft are located (refer to the map in Figure 4) at seven 
bases as follows (Table 10): 
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Table 10. Distribution of aircraft for example 


base 

A-10 

F-16 

F/A-18 

OEAH 

8 



KHARJ 


4 


OEJB 



4 

OEDF 


4 


OERY 

4 



OERK 


2 

2 

CV71 



8 


In this demonstration, we first run the static model with this information 
but consider only eight targets. Then we change the scenario and run the 
model again with the complete list of ten targets. This procedure simulates a 
situation in which new information becomes available moments after the first 
plan is disseminated. The added targets have more importance and must be 
attacked. Changes between the two solutions are highlighted in Table 11. 


Table 11. Results for two consecutives runs of the static model 
The first column shows the list of targets to be struck, but the first run only 
eight targets are available when the plan is generated. The second run 
simulates the situation where new information is suddenly available. Two targets 
with higher priority appear. The model generates a completely new plan._ 


target 

1st run (8 targets) 

2nd run (10 targets 
without persistence) 

i 

pckg 

base 

distance 

(NM) 

pk 

pckg 

base 

distance 

(NM) 

pk 

Safwan 

nOI 

OEAH 

353 

.95 

nOI 

OEAH 

353 

.95 

Al Asad 

nOI 

OEAH 

724 

.89 

nOI 

OERY 

678 

.89 

H2 

nOI 

OEAH 

773 

.95 

nOI 

OERY 

703 

.95 

H3 Airbase 

nOI 

OERY 

707 

.95 

NOT STRUCK 

H3 Highway 

nOI 

OERY 

719 

.95 

n05 

OERK 

705 

.91 

Wadi Al Khirr 

n04 

OERK 

495 

.90 

n04 

OERK 

495 

.90 

Tallil Air Base 

nOI 

OEAH 

442 

.95 

nOI 

OEAH 

442 

.95 

Rasheed 

n05 

OERK 

590 

.79 

NOT STRUCK 

Safwan2 

- 

- 

- 

- 

nOI 

OEAH 

353 

.89 

Safwan3 

- 

- 

- 

- 

nOI 

OEAH 

353 

.84 
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Figure 4. Location of targets and bases used in the examples 
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The result is that we have the important targets struck but five out of 
eight original allocations show major changes. We would prefer a new plan for 
which the addition of two targets has a lower impact. This is why we need a 
persistent model. 

The key outputs of the static model are optimal values of STRIKE(n,/J) 
which equals one if strike package n is assigned from site ; to target j. The 
persistent model will incorporate these results through the following 
parameters: 

prev(n,/,;') = previous optimal value of STRIKE(n,/,;) 

previj(/,j) = £prev(n,/,j) 

n 

prev_nogo(;) = previous optimal value of NOGO(;) 


A "persistence penalty” is applied to the new optimal value of 
STRIKE(n,/J) as follows 

route_pen persistence penalty for changing the base from which 
a target is attacked 

pkg_pen(n,n) persistence penalty for keeping the same base from 

which a target is attacked but changing the package 
used from n "to n 


pers_pen (n,i,j) = \ 


route _ pen if previj(/,;) = 0 and prev_nogo(;) = 0 
pkg _ pen(n, /T) if prev_nogo(;) = 0 and prev(/T,/, j) = 1 
0 otherwise 


The value of the pkg_pen parameter should be small when there is just a 
change in configuration, it should be larger if the number of aircraft of the 
same type increases, and should be yet larger if there is a change in aircraft 
type. 
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Finally, the objective function will incorporate: 

Minimize ^pers_pen (n,i,j) ■ STRIKE(n,/,;) 

n,i,j 


For the example, the objective function contains the following terms: 
Minimize ^tgtval(;) • NOGO(;') 

j 

+ ^ stress(n, i, j) ■ STRI KE(n, i, j) 

n,i,j 


+ Yj pers - pen(n, i, j) ■ STRI KE(n, /, j) 

n,i,j 


Table 12. Results for runs of the static model with persistence 
The same situation described in Table 11 is analyzed. But the model with 
persistence generates a plan with a reduced number of modifications of the 
original plan._ 


target 

1st run (8 targets) 

2nd run (10 targets 
with persistence) 

pckg 

base 

distance 

(NM) 

pk 

pckg 

base 

distance 

(NM) 

pk 

Safwan 

nOI 

OEAH 

353 

.95 

nOI 

OEAH 

353 

.95 

Al Asad 

nOI 

OEAH 

724 

.89 

nOI 

OEAH 

724 

.89 

H2 

nOI 

OEAH 

773 

.95 

nOI 

OEAH 

773 

.95 

H3 Airbase 

nOI 

OERY 

707 

.95 

NOT STRUCK 

H3 Highway 

nOI 

OERY 

719 

.95 

nOI 

OERY 

719 

.95 

Wadi Al Khirr 

n04 

OERK 

495 

.90 

n04 

OERK 

495 

.90 

Tallil Air Base 

nOI 

OEAH 

442 

.95 

nOI 

OEAH 

442 

.95 

Rasheed 

n05 

OERK 

590 

.79 

NOT STRUCK 

Safwan2 

- 

- 

- 

- 

nOI 

OERY 

380 

.89 

Safwan3 

- 

- 

- 

- 

nOI 

OERK 

363 

.84 


The solution based on the persistent model retains as much as possible 
from the original plan (Table 12). Only two modifications are necessary to 
accommodate the new targets. Table 13 shows that the persistent solution is 
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sub-optimal by only 1%; comparing Tables 11 and 12 shows that the persistent 
solution would be easier to implement. 


Table 13. Summary of the results for the static model example 



original 8 
targets run 

10 targets 
without 
persistence 

10 targets 
with 

persistence 

Major changes 

Targets not struck 

- 

2 

2 

Change of base 

- 

3 

0 

Total 

- 

5 

2 

Objective function value 

Value of targets not struck 

- 

180.00 

180.00 

Stress 

206.67 

151.75 

155.20 

Total 

206.67 

331.75 

335.20 


B. SUCCESSIVE EXECUTION OF THE DYNAMIC MODEL 

Before developing the persistent version of the dynamic model, we need 
to address re-execution of the model in successive time periods. 

Figure 5 represents the results of two successive executions. Time t 0 is 
the current time when the second problem is run. Time is the beginning of 
the first time period of the second run. If the first run plans a sortie for launch 
before and returns after t h then it must be regarded as a fixed decision in 
the second run. On the other hand, if the first execution plans a sortie that 
launches after ti, then this part of the plan can possibly be changed. The 
persistent model allows, but discourages, these changes. To account for 
previous plans, we define: 

prev_etd(n,/J) previous estimated time of departure of package n 
from site ; to attack target j, derived from the 
previous optimal values of STRIKE (n,i,j,w,nw,t) 

prev_flt(/,a,cj) flight time of aircraft a in configuration c from site 
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prev_nogo (j) 
fixed (j) 


; to target j in previous solution 

previous optimal value of the variable NOGO(;) 

equals one if it is too late to change plans for target 
j, zero otherwise 


This information and the updated time control parameters are used to 
specify fixed decisions and to re-derive the parameters dueback(/,a,t). 


present 
time (t 0 ) 


/resolution 


t-j r t2 t3 


Old ATO 


gap 


■ 4 —► 


horizon 


New ATO 


Figure 5. Super-imposition of a new plan onto an old plan 


This figure represents the results of two successive 
executions of the dynamic model. Horizontal lines represent the 
selected sorties. Sorties launched between t? and t 4 will be 
considered by both plans. 
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C. PERSISTENCE FOR THE DYNAMIC MODEL 

The persistent dynamic model can be re-applied many times, re¬ 
allocating strike packages considering new targets, threats and changes in 
priorities, weather or asset availability. 

The scenario for illustrating persistence in the dynamic model is the 
same as for the static model except that the packages are formed under the 
conditions described in Chapter III (refer to Table 7). The targets are given in 
Table 14, including the required time on target. 


Table 14. Target data example for the dynamic persistent model 


target 

latitude 

longitude 

target 

value 

time on 
target 

minimum 

damage 

deg 

min 

deg 

min 

Safwan 

30 

08 

47 

39 

100 

1000 

.80 

Al Asad 

33 

47 

42 

26 

150 

1015 

.75 

H2 

33 

21 

40 

35 

160 

1140 

.90 

H3 Airbase 

32 

55 

39 

44 

100 

1055 

.50 

H3 Highway 

32 

50 

39 

18 

180 

1100 

.40 

Wadi Al Khirr 

31 

25 

43 

11 

150 

1035 

.85 

Tallil Air Base 

30 

56 

46 

05 

200 

1200 

.85 

Rasheed 

33 

16 

44 

29 

80 

1040 

.50 


There are eight targets, it is 0800, and one hour is the minimum 
preparation time for the first departure. Figure 6 shows the result of the first 
run. 

Now suppose it is 0815 and we want to run the model again. All the data 
is the same, except for a change in the weather forecast from good to 
marginal, which affects the pk’s. Two of the launches from the 0800 solution 
were planned for 0900, so it is too late to change these decisions. These fixed 
decisions are indicated by asterisks in the 0815 solution, which appears in 
Figure 7. 
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Air Tasking Order (generated at: 800) 


TARGET 

TOT 

ETD 

A/C 

BASE 

DIST 

CONFIG 

WEAPON 

PK 

Safwan 

1000 

930 

2 

FA18 

OEDF 

194 

config02 

3 

MK-83 

0.81 

A1 Asad 

1015 

900 

2 

FA18 

CV71 

587 

configOl 

2 

MK-84 

0.75 

H2 

1140 

900 

2 

A10 

OERY 

703 

configOl 

4 

MK-82 

0.90 

H3 Airbase 

1055 

915 

1 

FA18 

CV71 

679 

configOl 

2 

MK-84 

0.50 

H3 Highway 

1100 

930 

1 

FI 6 

OERK 

705 

config02 

2 

MK-84 

0.45 

Wadi A1 Khirr 

1035 

915 

2 

FA18 

CV71 

453 

config02 

4 

MK-83 

0.88 

Tallil Air Base 

1200 

1030 

2 

A10 

OERY 

432 

config02 

2 

MK-83 

0.85 

Rasheed 

1040 

930 

1 

FA18 

CV71 

479 

configOl 

2 

MK-84 

0.50 


Figure 6. Solution to the dynamic model run at 0800 

The output for the dynamic model is an ATO with the 
following information: name of the target, required time on 
target (TOT), estimated time of departure (ETD), number and 
type of aircraft participating on the attack (A/C), distance (in 
nautical miles) from the base to the target, configuration to be 
used, number and type of the weapon to be delivered and the 
expected probability of kill. 


Air Tasking Order (generated at: 815) 


TARGET 

TOT 

ETD 

A/C 

BASE 

DIST 

CONFIG 

WEAPON 

PK 

Safwan 

1000 

930 

2 

FA18 

OEDF 

194 

config02 

3 

MK-83 

0.81 

A1 Asad 

1015 

900 

2 

FA18 

CV71 

587 

configO1 

2 

MK-84 


H2 

1140 

900 

2 

A10 

OERY 

703 

configO1 

4 

MK-82 


H3 Airbase 

1055 

915 

1 

F16 

OERK 

693 

configOl 

2 

MK-83 

0.53 

H3 Highway 

1100 

930 

1 

FA18 

CV71 

699 

config02 

2 

MK-83 

0.46 

Wadi A1 Khirr 

1035 

930 

2 

FA18 

CV71 

453 

config02 

4 

MK-83 

0.88 

Tallil Air Base 

1200 

1030 

2 

A10 

OERY 

432 

config02 

2 

MK-83 

0.85 

Rasheed 

1040 

915 

1 

F16 

OEDF 

567 

configOl 

2 

MK-83 

0.53 


Figure 7. Solution to the dynamic model run at 0815 without persistence 

Sorties marked with an asterisk are fixed. Differences from 
the 0800 solution are highlighted in bold. 


The pk’s for the fixed sorties are left blank. This is because the 
previously assigned weapon may be ineffective in the marginal weather state, 

but it is too late for this run of the model to consider changing weapons. 
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The version of the dynamic model used to obtain this solution does not 
encourage persistence. Different bases and aircraft are assigned for three of 
the six targets not constrained by the earlier run. This is a major disruption of 
the plan in only 15 minutes. 

The next step is to implement the persistence penalties in the same way 
we did for the static model in section A. Penalty parameters route_pen and 
pkg_pen (n,n't) are applied to the new optimal value of STRIKE (n,i,j,w,nw,t), 
for changing routes and packages. We define the parameter pkg_pen (n,n't) 
using the following criteria. Assume package n uses nac aircraft of type a in 
configuration c; similarly define n' a', nac', and c'. Then pkg_pen(n,n jt) = 0 if 
n = n' otherwise, and pkg_pen (n,n't) is the sum of penalties as follows: 

pi persistence penalty for different aircraft type 

p2 persistence penalty for smaller number of aircraft 

p3 persistence penalty for greater number of aircraft 

p4 persistence penalty for different configuration 

prev_forecast(t) previous weather forecast for period t 


pi if a* a' 

p2 ■ (nac' - nac) if a * a' and nac < nac' 
p3 ■ (nac - nac') if a * a' and nac > nac' 
p4 if a*a',c*c' and forecast (t) prev_ forecast (t) 


In other words, there is no penalty when the same package is assigned in 
the second solution as in the first but there are cumulative penalties for using 
different types of aircraft, different number of aircraft or different weapon 
configurations. However, there is no penalty for changing weapons when the 
weather changes. 
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The persistence penalties applied to STRIKE (n,i,j,w,nw,t) variables in 
the dynamic model are exactly the same as we had for the static model: 


pers_pen (n,i,j,t) = 

route _ pen if previj(/,;) = 0 and prev_nogo(;') = 0 
<pkg_pen (n,n',t) if prev_nogo(;) = 0 and prev(n',/,;) = 1 
0 otherwise 


Finally, the dynamic persistent model has the following objective 
function: 


Minimize 


^ pers_pen (n,i,j,t) ■ STRIKE (n,i,j,w,nw,t) 

n>i>j>w,nw,t 

+£tval(j) NOGO(j) 

1 


+ I 


dpen(/, a, j) + pkpen(n, j, w, nw) 
+attrition(n, j) + dnac(n, a) ■ acval(a) 

STRIKE(n,/,;,w,nw, t) 


Using this model, we are able to create the persistent solution presented 
in Figure 8. 


Air Tasking Order (generated at: 

815)/p 






TARGET 

TOT 

ETD 

A/C 

BASE 

DIST 

CONFIG 

WEAPON 

PK FIX 

Safwan 

1000 

930 

2 

FA18 

OEDF 

194 

config02 

3 

MK-83 

0.81 

A1 Asad 

1015 

900 

2 

FA18 

CV71 

587 

configO1 

2 

MK-84 

• k 

H2 

1140 

900 

2 

A10 

OERY 

703 

configO1 

4 

MK-82 

k 

H3 Airbase 

1055 

915 

1 

FA18 

CV71 

679 

config02 

2 

MK-83 

0.57 

H3 Highway 

1100 

930 

1 

FI 6 

OERK 

705 

configOl 

2 

MK-83 

0.53 

Wadi A1 Khirr 

1035 

915 

2 

FA18 

CV71 

453 

config02 

4 

MK-83 

0.88 

Tallil Air Base 

1200 

1030 

2 

A10 

OERY 

432 

config02 

2 

MK-83 

0.85 

Rasheed 

1040 

930 

1 

FA18 

CV71 

479 

config02 

2 

MK-83 

0.57 


Figure 8. Solution to the dynamic model run at 0815 with persistence 


In this persistent solution, only three minor changes occur; 
they are noted in bold type. 
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All bases remain the same for this solution. The three configuration 
changes were expected because in the original plan the MK-84 weapon is used 
many times but is less effective under the new weather forecast. As before, all 
targets are struck. Table 15 summarizes and compares the results for each run. 


Table 15. Summary of the results for the dynamic model example 



original run 

subsequent 
run without 
persistence 

subsequent 
run with 
persistence 

Major changes 

- 

6 

0 

Change of base 

- 

3 

0 

Change of aircraft type 

- 

3 

0 

Objective function value 

33.978 

34.271 

35.726 


The objective function values reported in Table 15 exclude persistence 
penalties and the terms attributed to fixed decisions. The original run has the 
best objective function value, because it assumes better weather. 
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V. CONCLUSION 


This research has developed models for optimally composing air strike 
packages and assigning these packages to strike a set of prioritized targets. A 
package consists of aircraft from various bases or carriers along with 
appropriate weapons for the assigned target. The models ensure that the 
aircraft have sufficient range, time on target and the right weapons so that a 
sufficiently high probability of kill is achieved for each assigned target. 

Two basic models are created, a static one and a dynamic one. The 
static model covers a short time frame during which an aircraft would fly at 
most one sortie. This model allows the creation of heterogeneous strike 
packages with different sorts of strike aircraft. The dynamic model covers a 
longer time frame during which an aircraft may be involved in multiple strikes. 
For computational reasons, this model only considers homogeneous packages, 
i.e., packages containing a single type of aircraft. 

Both models have "persistent” variants. A persistent model is important 
when a strike plan is already in place, new high-priority targets arise, and a 
new plan must be developed. The persistent model creates a new strike plan 
that disrupts the old one as little as possible - this is important because much 
time-consuming work has probably already been invested in implementing the 
original plan - yet is still near-optimal in the mathematical sense. 

All models precompute nonlinear probabilities of kill for potential strike 
packages. Binary variables represent the assignment of these potential 
packages to targets so that a mixed-integer linear program results. An instance 
of the dynamic model with three weapon types, 36 aircraft from seven bases 
and 100 potential targets is solved in less than three minutes on a personal 
computer. 
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A. REFINEMENTS 

1. Weapons Effectiveness Calculations 

Our models use simple, fast and well-known algorithms to calculate the 
probability of kill for a combination of assets, from information such as a 
weapon’s lethal radius and circular error probable. However our algorithms do 
not cover all weapon and target types, and depending on the nature of the 
application, other algorithms should be used. The Joint Technical Coordinating 
Group for Munitions Effectiveness (JTCG/ME) Program, for example, develops 
analytical methods for measuring and predicting munitions effectiveness [USAF 
1998]. This group has also produced a large body of scientifically valid data 
related to specific weapons, munitions, and appropriate targets. Their models 
are more sophisticated than ours and include such details as aircraft 
capabilities and configurations; target characteristics, such as size, shape, and 
hardness; and delivery parameters such as altitudes, speeds, and dive angles. 
These models and algorithms, or other options, should be considered for 
refining our models’ probability-of-kill calculations. 

2. Weather Effects 

The way we handle weather is perhaps too simple. We only modify 
weapons-effectiveness calculations based on three weather states, but the real 
situation is more complex. Even slight changes in humidity, ambient light, or 
intensity of precipitation can significantly degrade weapons systems 
performance. Sometimes, a pilot will detect a target, but his weapons systems 
cannot "see” the target because of weather conditions. The results are an 
increase in ammunition expenditure, and degradation in time-on-station 
performance and mission accomplishment [Cohen 1993]. Consequently, we 
believe that more refined weather information could be usefully incorporated 
into our models. 
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B. FURTHER RESEARCH PROJECTS 

There are several potential projects that can be pursued to extend the 
models presented in this thesis. This thesis has only addressed the strike 
aircraft of an air strike package. Additional work, in both the static and 
dynamic models, will be required to include specialized escort aircraft and 
aircraft for suppression of enemy air defenses. However, our basic model 
paradigm should extend directly: we precompute the set of potential strike 
packages and the model optimally selects a feasible subset of those packages 
and assigns them to targets. 

An important extension of the dynamic model would allow the model to 
choose the optimal time on target. The target list should specify a time window 
for the strike and let the model choose time and the best assets to use. A more 
difficult extension of the dynamic model would incorporate heterogeneous 
packages. 

Another challenging project would permit the redirection of previously 
committed sorties. In an extension of the current dynamic model, aircraft 
enroute to targets could possibly be redirected to a just-identified, higher- 
priority target, or to a secondary target if the original target is no longer 
appropriate for some reason (because of changes in local weather conditions, 
new battle-damage information shows the target has already been destroyed, 
etc.). New variables will be necessary to monitor each aircraft’s position and to 
determine which targets are close enough to be struck and how much air 
refueling is needed. 

Another extension of this work would allow strike package sorties to 
strike multiple targets. One way this could happen is if the package stays intact 
and visits two targets with weapon loadouts sufficient for both. A more 
difficult scenario to model and solve would allow a package to strike a primary 
target and then divide into "sub-packages” to strike secondary targets. 
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APPENDIX A. CALCULATIONS OF THE PROBABILITY OF KILL 


This appendix presents the formulation used to calculate the probability 
of kill of the packages created for the static and dynamic models. The 
reference for this subject is "Notes on Firing Theory” [Washburn 2002]. 

To calculate the probability of kill we should consider solving an 
independent sequence of n aircraft delivering sets of weapons arranged in 
dependent shots. The independent part is easily formulated as 

Pm of n indep saivos = 1 - [(I ~ ^) • 0 - ^)• • • 0 ~ )] = 1 - fl 0 - ^) where P kj is 

7=1 

the probability of kill of the 7-th aircraft. If the package has only one type of 
aircraft and each aircraft delivers the same quantity and type of weapons, the 
formula simplifies to 

Pm of n identical in d e P saivos = 1 - 0 - P S where n is the number of aircraft and 
P k is constant among them. 

In both cases we still need to calculate the probability of kill for each 
aircraft. In the case where we have different types of aircraft and weapons in 
the package, the problem is more complex: each aircraft has its own delivery 
profile (angle, velocity, etc.) for some type of weapon that will be different 
from another aircraft, even if all aircraft are launching the same type of 
weapon. Those particularities will be reflected as different parameters for 
each aircraft-weapon pair. 

However, we can have all the values of pk calculated in advance for 
each package against each type of target. For the purpose of this thesis, we 
apply a simplified assumption described in Washburn [2002]. We consider the 
weapon having a "diffuse Gaussian” (DG) damage function which, more 
realistically, should be applied only to weapons that kill by fragmentation. The 
DG assumption produces a simple expression for P kj : 
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[ s I 

m-p, (U,V)) is the probability of kill of a collection of s 

v <j =1 J 

shots from a given aircraft, where the expected value is with respect to the 
normal distribution of the random common error ( U,V ). Each aircraft’s shot can 
potentially have distinct parameters for dispersion and lethality which allows 
us to determine the result for a heterogeneous package. 

For the homogeneous package, we start with the simple case of 
calculating the probability of kill of a single shot: 

P m = 1 - (.5) (R2/CEp2) where R is the lethal radius of a "cookie cutter” air- 
to-ground weapon and CEP is the circular error probable (CEP = 1.1774 o 
where o is the standard deviation associated to the weapon). 

The first parameter is dependent on the type of weapon and the nature 
of the target. The second is a function of the weapon’s aerodynamics and 
aircraft’s delivery profile. 

Now, b weapons launched from the same platform represent dependent 
shots subject to the same bias and pattern of the distribution of the hits on the 
ground, so we use the result that Washburn calls the "confetti approximation:” 

hR 2 

P k = 1 - (1 + VTT) ■ exp(-VTT) where Z = 

2(7 
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